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Investigator: Dr. Stephen Walch 


A manuscript describing the calculations on the 1 CH 2 + H 2 0, H 2 + HCOH, 
and H 2 + H 2 CO product channels in the CH 3 + OH reaction, which were de- 
scribed in the last progress report, has been accepted for publication in J. Chem. 
Phys., and a copy of the manuscript is included in the appendix. The production 
of 1 CH 2 in this reaction is important in hydrocarbon combustion since 1 CH 2 is 
highly reactive and would be expected to insert into N 2 , possibly leading to a new 
source for prompt N0 X (vide infra). 

The appendix also contains a copy of a paper (J. Chem. Phys., in press), which 
discusses the barriers for 

H + CH 2 0 <-> ClhO <-> CH 2 OH 

These species are also important in hydrocarbon combustion and the computed 
barrier heights are consistent with experimental work (F. Temps et al.) in con- 
trast to earlier studies which obtained barrier heights which are too large. 

During the last six months new calculations have been carried out for the NH 2 
+ NO system, which is important in the thermal de-NO x process. Fig. 1 shows 
a schematic of the computed potential energy surface (PES) for this system. Ad- 
dition of NH 2 to NO leads to initial formation of an NH 2 NO species (1). This 
species is found to be nonplanar, with a pyramidal geometry about the NH 2 . For- 
mation of N 2 + H 2 0 proceeds via a complex series of steps starting with a 1,3- 
hydrogen shift to give a planar HNNOH species (2), in which the substituents are 
trans with respect to the NN bond and cis with respect to the NO bond. Rota- 
tion about the NO bond leading to (3) involves only a small barrier. However, 
rotation about the NN bond to give (4), which is the favorable configuration for 
elimination of H 2 0, requires breaking an NN 7r bond and this pathway involves 
a very high barrier. Instead conversion of (3) and (4) is found to proceed via a 
planar saddle point which may be described as an inversion process. Subsequent 
elimination of H 2 0 occurs via a planar saddle point with an approximately rect- 
angular arrangement of H-N-N-O. (See Fig. 2) The other product channel is HN 2 
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+ OH, which is expected to occur from (2), (3), or (4) with no barrier other than 
the exothermicity. Since the HN 2 + OH asymptote is at about the same energy 
as the bottle neck for formation of H 2 0 + N 2 , the calculations suggest that both 
product channels are possible. An accurate estimate of the product branching ra- 
tio requires dynamics calculations ( vide infra). 







N i-t di - 0 


Fig. 1. Schematic of the PES for NH 2 + NO. 

In a collaboration with Schatz (Northwestern) studies are being undertaken to 
determine the appropriate dynamics methods for determining product branching 
ratios in these complex reactions. As initial systems for study we have selected 
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Fig. 2. Saddle point geometry for HNNOH (4) -+ H 2 0 + N 2 . 

the HNO and HN + NO systems. In the first case the reaction to be studied is 
+ o —>• N 4- OH or H + NO. Characterization of this system requires global 
PES’s for the 1 A', 3 A", and 1 A" surfaces of H + NO. Current calculations on this 
system involve 19 values of the Jacobi angle, up to 10 values of R, and 4-6 values 
of r, or about 800 points to characterize the HNO and HON regions for three sur- 
faces. Currently about 2/3 of the points have been completed. The global PES 
will be based on spline fits of the Morse parameters in r as was done previously 
for HN 2 . It wih also be necessary to represent the N + OH and NH + O chan- 
nels, however, it is expected that this region of the surface can be represented 
with a simple functional form which requires few additional points. Once the sur- 
face is complete it will be possible to test various dynamical approximations. 

For the NH + NO reaction, experiments have focused on H + NNO which leads 
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to NH + NO and N 2 + OH as major product channels. It is planned to spice to- 
gether this surface from the known triatomic potentials (HNN, HNO, and NNO) 
plus information on the four atom system which is already available from work 
done previously under this grant ( see the previous progress report). It is antic- 
ipated that a grid of points in the N 2 + OH region will also be needed to help 
characterize the portion of the surface involved in energy release. These studies 
should help to understand interesting experimental results (Wittig et al.) which 
show substantial internal energy ( presumably vibrational energy) in the N 2 prod- 
uct. 

In addition to continuation of the work discussed above future work will in- 
clude studies of the 1 CH 2 + N 2 reaction mentioned above. The initial product 
on the singlet surface is expected to be a cyclic CH 2 N 2 isomer (C 2xr symmetry). 
This species should have a fairly small bond energy ( probably of the order of 
magnitude of 20 kcal/mol.). The cyclic isomer could open to form the singlet 
ground state of diazomethane or via a singlet-triplet surface crossing could lead 
to CH2N(doublet) + N( 4 S). The latter process is expected to be endothermic by 
analogy to the isoelectronic reaction, 

OCD) + N 2 - NO + N{ 4 S) (1) 

which is endothermic by a 30 kcal/mol. A more likely low barrier process is hy- 
drogen abstraction from the initially formed cyclic CH 2 N 2 to give a cyclic CHN 2 
species (C 2 „ symmetry). This is the same species proposed by Yarkony in the 
analagous reaction 

CH( 2 U) + N 2 -» HCN + N( 4 S) (2) 

This process requires a singlet-triplet surface crossing, however, Yarkony’s calcu- 
lations show that the intersection of the two surfaces occurs at approximately the 
geometry of the cyclic CHN 2 species. In either case, the N atoms can enter into 
the second step of the thermal Zeldovich mechanism i.e. 

N + 0 2 ^N0 + 0 (2) 
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which leads to N0 E production. Initial work on this system will focus on the en- 
ergetics of the species involved (i.e cyclic CH 2 N 2 , diazomethane, and CH 2 N + 
N( 4 S) ). CH 2 N if formed is expected to be highly reactive to 0 atoms formed 
by reaction (2). The expected products are 1 CH 2 + NO. Thus, if production of 
CH 2 N is an important pathway the 1 CH 2 species serves a catalytic role in the net 
process 

N 2 + 0 2 - 2 NO (3) 

Another system which will be studied is H + HNO. This reaction may be im- 
portant in experiments to simulate supersonic combustion where heated air is 
used resulting in formation of NO via the thermal (Zeldovich) process. Here HNO 
may be an intermediate in the recombination reaction 

H + H + NO -» H 2 + NO 

This system ties in with previous work on the NH 2 0 system and with current 
work on HNO. 
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Investigator: Dr. Eugene Levin 


During the time period from 1 July 1992 to 31 December 1992, the determina- 
tion of the transport properties of gases (e.g., viscosity, diffusion, thermal conduc- 
tivity, etc.) from the collision cross-sections of the constituent species was com- 
pleted for the mixtures H-H 2 , H-H, H-Ar, and H-H 2 0. This work was carried out 
in collaboration with NASA scientists St allcop, Partridge (and others) and four 
papers documenting this work are in preparation for submission to the Journal of 
Chemical Physics. 

The results for H-H 2 O were particularly significant in three respects. Since the 
water molecule is in a C2V configuration rather than linear, a specific plane is 
determined. Hence, it was necessary to determine representative cross-sections 
for a very large number of collision approach directions between the hydrogen 
atom and the water molecule representing both in-plane angular variations as 
well as out-of-plane relative orientations. The only previous attempt to apply our 
method to other than binary collisions was a study of H-N 2 interactions (“H-N 2 
interaction energies, transport cross sections, and collision integrals”, by Stallcop, 
Partridge, Walch and Levin, JCP 97 (5), September 1992, pp. 3431-3436). That 
study did involve consideration of the various angles of approach, but since the 
N 2 molecule is linear, there was one less degree of freedom than in the case of II- 
H 2 0. For water, over 50 different collision geometries were considered in order to 
obtain a reliable weighted mean cross section. 

A second important aspect of the H-II 2 0 study was the discovery that the av- 
erage cross-sections and collision integrals resulting from the weighted mean of 
the individual angular collision orientations differed very little from the results 
obtained by using an average potential energy representation. A similar result 
was found in the H-N 2 study cited above but it was unexpected for the more 
complex case of H-H 2 0. If this proves to be a general result, it would simplify 
greatly the determination of cross sections and transport properties for inter- 
actions of atoms or ions with larger molecules. It may also apply for molecule- 
molecule collisions. 
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The third important feature of the H-H 2 O study addressed the question of 
whether the semi-classical approximation used to calculate the phase shifts re- 
tained sufficient accuracy for small values of the reduced mass, /i, particularly at 
low energy levels. For H-H 2 0 the reduced mass is approximately 0.947 and en- 
ergy levels as small as E=0.0001 hartrees were used in the study. The error in 
cross sections due to the semi-classical approximation was determined by direct 
integration of the radial Schroedinger equation to determine the phase shifts for 
a range of energy values using the average potential energy representation for 
the interaction of H with H 2 0. It was found that the error in the cross sections 
diminished from about 12% at E=0.0001 to less than one percent at E=.002. 

The resulting error in the collision integrals decreased from about 4% at T=100 
(Kelvin) to less than 1% at T=400 K. It was concluded that the semi-classical 
method remained valid except for extremely low energies and the correspondingly 
low temperatures. 

As a result of the above findings for H-H 2 0 regarding the range of validity of 
the semi-classical method, all of our current work in the determination of cross 
sections and collision integrals now uses a quantum mechanical calculation for 
low energies and small reduced masses fi and only uses the semi-classical approx- 
imation above the energy level for which the difference is less than a fraction of 
1%. Our earlier work on the transport properties of the constituents of air was 
rechecked to see if the limitations of the semi-classical method had resulted in er- 
roneous results, however it was found that the published results were not affected 
since the reduced masses for interactions of 0-0, N-N and N-0 (and their ions) 
were much larger than collisions involving hydrogen. In all of these cases the dif- 
ferences between the semi-classical and quantum mechanical results were negligi- 
ble. 

Since the direct integration of the radial Schroedinger equation to obtain the 
quantum mechanically determined phase shifts is very time-consuming, a much 
more efficient method was developed by Schwenke which combines the techniques 
of R-matrix propagation and Richardson extrapolation. The algorithm is des- 
ignated the HEMPR method (h-extrapolated Magnus propagator R-matrix). 

The details of this new algorithm will be contained in a joint paper by Levin, 
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Schwenke and Stallcop which compares the results of quantum mechanical and 
semi-classical determinations of cross sections and collision integrals for various 
energy levels and values of the reduced mass fi. 

Since hydrogen is at least one of the collision partners, the reduced masses are 
small and consequently the results obtained for H-H, H-H2, and H-Ar all were 
based on the HEMPR method at low energies and the semi-classical method af- 
ter the differences became negligible. The results were found to be in quite good 
agreement with available highly accurate (but limited energy range) experiments. 
The extension of the cross section determination over a broad spectrum of en- 
ergy levels permitted the transport properties to be presented over a wide range 
of temperatures of interest to NASA. 

In addition to completing the documentation and publication of the work de- 
scribed above, the major goals of next time period are to undertake the study of 
molecule-molecule collisions (specifically, H2-H2, H2-N2, and N2-N2) and to com- 
plete an Air Reference Model. The molecule- molecule study is a significant ad- 
vance over our prior atom-atom and atom-molecule work. Of particular interest 
will be to determine whether the average cross sections obtained by combining 
the cross sections for all the various possible collision orientations is or is not the 
same as the resulting cross sections using a mean potential. The completion of 
the Air Reference Model is awaiting the determination of the contributions from 
the collisions N-e _ and 0 -e - . 
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Theoretical Characterization of the Reaction 

CH 3 + OH — > CH 3 OH — ► products : 

The 1 CH 2 + H 2 0, Ho + HCOH, and H 2 + HoCO Channels 

Stephen P. Walch a 
ELORET Institute 
Palo Alto, Ca. 94303 

Abstract. The potential energy surface (PES) for the CH3OH system has been char- 
acterized for the : CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO product channels 
using complete-active-space self-consistent-field (CASSCF) gradient calculations to 
determine the stationary point geometries and frequencies followed by CASSCF/ 
internally contracted configuration-interaction (CCI) calculations to refine the ener- 
getics. The 1 CH 2 + H 2 0 channel is found to have no barrier. The long range inter- 
action is dominated by the dipole-dipole term, which orients the respective dipole 
moments parallel to each other but pointing in opposite directions. At shorter sep- 
arations there is a dative bond structure in which a water lone pair donates into the 
empty a ,r orbital of CH 2 . Subsequent insertion of CH 2 into an OH bond of water 
involves a non-least-motion pathway. The H 2 + HCOH, and H 2 + H 2 CO pathways 
have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect 
to CH 3 + OH. From comparison of the computed energetics of the reactants and 
products to known thermochemical data it is estimated that the computed PES is 
accurate to =t 2 kcal/mol. 

a Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 


The CH 3 + OH reaction has at least six possible product channels: 


CHz +OH + M 

-> CH 3 OH + M 

(la) 

CH 3 4- OH — > 

CH 2 OH + H 

(16) 

CH 3 + OH — > CH z O + H 

(lc) 

CH Z + OH —* 

1 CH 2 + H 2 0 

(Id) 

CH 3 + OH — * 

HCOH + H 2 

(lc) 

CH Z + OH - 

* H 2 CO + h 2 

(1/) 


The role of reaction Id has been controversial. The room temperature rate 
for the reverse of reaction Id has been measured by Hatakeyama et al. [1] 
as ~ 3 x 10 -12 cm 3 molecule -1 cm - 1 , while Hack et al. [2] obtained 
3.5 x 10 -n cm 3 molecule -1 cm -1 , or about an order of magnitude faster. The latter 
rate is approximately gas kinetic and implies no barrier. Using currently accepted 
heats of formation [3], and a singlet-triplet splitting for methylene of 9 kcal/mol, 
reaction Id is exothermic by 0.7 kcal/mol at OK. Thus, the forward reaction is also 
expected to be very fast. 

Dean and Westmoreland [4] have used a variant of Rice-Ramsperger-Kassel- 
Marcus (RRKM) theory called QRRK theory to model the product distributions 
in the CH 3 + OH reaction. In this work the parameters for lb and lc were based 
on estimated rates for the reverse reaction. The parameters for reaction If were 
taken from calculations [5] and the rate for reaction Id was taken from Ref. 1. 
According to these studies, at room temperature and moderate pressure CH 3 OH is 
the dominant product, while at flame temperatures reaction lb is thought to take 
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over. The HCOH + H2 channel does not appear to have been considered in these 
studies, though theory [5] indicates essentially no barrier with respect to CH 3 + 
OH. This study indicated that production of 1 CH 2 (reaction Id) is a minor channel. 
By contrast a model proposed by Pilling and coworkers [6] which makes use of the 
rate for the reverse of reaction Id due to Hack et al. [2] indicates that reaction Id 
is the dominant channel above room temperature. 

Recently Smith [7] has also reported RRKM calculations for CH3 + OH. These 
calculations as well as the work of Pilling et al. [6] have indicated a need for more 
accurate potential energy surface (PES) information for the 1 CH 2 + H 2 0, H 2 + 
HCOH, and H 2 + H 2 CO product channels. The most accurate previous theoretical 
study of these channels in the CH3OH system was carried out by Harding, Schlegal, 
Krishnan, and Pople [5] using M0ller-Plesset perturbation theory with a 6-311G** 
basis set. Although these calculations were carefully carried out, there is probably 
considerable uncertainty in the energetics; by current standards, both the basis set 
and treatment of electron correlation can be improved upon. More recently the bond 
dissociation energies of CH 3 OH have been computed by Bauschlicher, Langhoff, and 
Walch [8] using the modified coupled-pair functional method. Similar calculations 
were carried out by Pople and co-workers [9] using the G2 method, which includes 
some empirical corrections. These calculations accurately determined the heats of 
formation of the CH 3 0 and CH 2 OH species, but they did not consider the portions 
of the PES leading to the 1 CH 2 + H 2 0, H 2 + HCOH, and H 2 + H 2 CO product 
channels. Thus, these regions of the PES are reexamined here. 

Qualitative features of the potential energy surfaces are discussed in Sec. II, the 
computational method is discussed in Sec. Ill, the results are presented in Sec. IV, 
and the conclusions are given in Sec. V. 
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II. Qualitative Features. 

CH 2 has two low-lying states. The 3 B a ground state which will be drawn as 



has two orthogonal high-spin coupled singly occupied orbitals, while the state 
which is « 9 kcal/mol higher, is drawn as: 


H 

H 



i.e., the two orbitals corresponding to the C lone pair are singlet paired. If these 
orbitals are solved for self consistently in a generalized valence bond wavefunctaon 
[10] the overlap integral between them is » 0.7. Thus, the singlet state of CH 2 may 
be characterized as a singlet biradical, but the orbitals of the lone pair have a sub- 
stantial overlap, which must be maintained while inserting into a single bond if the 
process is to occur without a substantial barrier. In terms of multi configurational 
self-consistent-field (MCSCF) theory, the biradical character in CH 2 arises because 
of a near degeneracy effect between an sp hybrid lone pair and an empty 2p-like a" 

orbital. 

The insertion of 1 CH 2 into a bond pair occurs via a non-least-motion pathway. 
This process is completely analagous to the non-least-motion addition of the 2 II 
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state of CH to H 2 , as discussed by Dunning and Harding [10]. In that case the CH 
and H 2 approach each other with the bond axes parallel and for this orientation 
the C lone pair and H 2 bond pair orbitals are able to evolve into two CH bonding 
orbitals without breaking either bond. These orbital changes are consistent with 
the orbital phase continuity principle arguments made by Goddard [11]. 

In the case of 1 CH 2 + H 2 0 the long range interaction is dominated by the dipole- 
dipole term, which leads to an initial approach with the dipoles parallel to each other 
but pointed in the opposite direction. At shorter separation there is a dative bonded 
structure where one lone pair on H 2 O donates into the empty C 2p-like orbital of 
CH 2 . This structure is a minimum on the PES. Insertion of CH 2 into one of the 
H 2 O bonds requires rotating H 2 O such that one OH bond is in the plane defined by 
O and the bisector of L HCH. This orientation is similar to the orientation favored 
for non-least-motion insertion of CH into H 2 . 

Insertion of 1 CH 2 into H 2 O involves very little barrier. The insertion of hydroxy 
methylene (HCOH) into H 2 , however, is found to exhibit a barrier of 13.9 kcal/mol. 
The lower reactivity of HCOH as compared to CH 2 results from delocalization of 
an 0 lone pair into the empty C 2p orbital. The resultant exclusion effect reduces 
the 2s — + 2p near-degeneracy effect and thus results in less biradical character in 
the substituted carbene. 

III. Computational Details. 

Two different basis sets were used in this work. For the CASSCF gradient calcu- 
lations the polarized double-zeta set of Dunning and Hay [12] was used. The basis 
set for C and 0 is a (9s5p)/[3s2p] basis augmented by a single set of 3d functions 
with exponents of 0.75 and 0.85 for C and O, respectively. The H basis is (4s) / [2s] 
augmented with a single set of 2p functions with exponent 1.00. The basis set 
used in the Cl calculations is the Dunning correlation consistent triple-zeta double- 
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polarization basis set [13]. This basis is [4s3p2dlf] for C and O and [3s2pld] for H 
and is described in detail in Ref. 13. 

In the CASSCF calculations the electrons in the bonds which are being made 
or broken were included in the active space. In the case of the 1 CH2 + H2O and 
H2CO + H2 channels, there were six active electrons and six active orbitals. The 
correlated electrons correspond to the CO bond, the OH bond, and one CH bond 
of CH3OH and evolve to the C lone pair of CH 2 and the two OH bond pairs in the 
case of 1 CH 2 + H 2 0 and the CO a and 7 r bonds and the H 2 bond in the case of 
H 2 CO + H 2 . The remaining 12 electrons are inactive in the CASSCF calculation. 
In the case of the HCOH + H 2 channel, there were four active electrons and four 
active orbitals. The correlated electrons correspond to two CH bonds in CH3OH 
and to the C lone pair in HCOH and the H 2 bond pair in the HCOH + H 2 limit. 
The remaining 14 electrons were inactive in the CASSCF calculation. In generating 
the set of reference configurations, no more than two electrons were permitted in 
the weakly occupied CASSCF orbitals. All but the 0 Is and C Is electrons were 
correlated in the CCI calculations. 

The CASSCF/gradient calculations used the SIRIUS/ABACUS system of pro- 
grams [14], while the CCI calculations were carried out with MOLPRO [15,16]. 
Most of the calculations were carried out on the NASA Ames Cray Y-MP; although 
some of the CCI calculations were carried out on the NAS facility Y-MP. 

IV. Discussion. 

Table I shows computed energetics for the portions of the CH3OH surface which 
were considered in this work. The computed energetics were obtained from the CCI 
energies (including the multireference analogue of the Davidson correction [17]) plus 
the zero-point energies from the CASSCF calculations. (See Tables Ila-IIc.) Thus, 
these energetics should be compared to experimental values at 0 K. The computed 
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frequencies and rotational constants for the saddle points for CH2O + H2, HCOH 
+ H2, and 1 CH 2 + H2O, ( denoted as CH2O-H2, HCOH-H2, and CH2-H2O, re- 
spectively) and for the CH2 + H2O dative bonded structure (denoted as CH2.H2O) 
are also given in Table III. The stationary point corresponding to CH2.H2O is a 
minimum on the PES, but there is one very small frequency ( 26 cm *), which 
corresponds to a hindered rotation of the CH2 and H2O with respect to each other. 

The computed energy separations discussed here, in each case, involve breaking 
two bonds and forming two new bonds; thus, the errors in the individual bond 
strengths cancel and the computed energetics are expected to be accurate. However, 
in Ref. 2 it was shown that, for calculations of about the same quality as reported 
here, the error in the C -0 bond strength in CH3OH is 6.5 kcal/mol. Thus, in order 
to compute energies with respect to CH 3 + OH, the experimental 0 K value of 
90.2 kcal/mol [ 18 ] was used for the C -0 bond strength. The locations of the H + 
CH 3 0 and H + CHoOH asymptotes were taken as the best-estimate values from Ref. 
2 . This places H + CH 3 0 and H + CH 2 OH at 14.8 kcal/mol and 6.0 kcal/mol above 
CH 3 + OH, respectively. The experimental locations of a CH2 + H2O and CH2O + 
H2 with respect to CH 3 OH were derived from the JANAF [ 3 ] heats of formation of 
H 2 , CH 2 0 , CH 3 , OH, H 2 0 (g), and 3 CH 2 , plus a singlet-triplet splitting in CH 2 of 

9.0 kcal/mol and the value for the C -0 bond strength in CH3OH given above. 
From Table I it is seen that the computed CH 3 OH — * 1 CH2 + H2O separation 

is 0.7 kcal/mol smaller than experiment, while the computed CH3OH — ► CH2O + 
H2 separation is 0.6 kcal/mol smaller than experiment. It is also seen that the 
computed results of Harding et al. [ 5 ] are 5.4 kcal/mol larger and 2.3 kcal/mol 
smaller, respectively, for the same separations. The computed barrier heights for 
CH3OH — ► HCOH + H 2 and CH 3 OH -* CH 2 0 + H 2 obtained in Ref. 5 are 

6.0 kcal/mol and 4.6 kcal/mol larger than the barrier heights obtained in the present 
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work. These differences presumably reflect a combination of the larger basis set 
and more extensive correlation treatment in the present calculations. Based on the 
comparison to experiment for the 1 CH2 + H2O and CH2O + H2 asymptotes it is 
reasonable to assign error bars of ~ ± 2 kcal/mol to the computed energetics in 
this work. 

The computed energetics for the CH3OH system are also shown schematically 
in Fig. 1. From Fig. 1 it is seen that only the CH 2 0 + H 2 channel exhibits a 
barrier ( 1.7 kcal/mol.) with respect to CH3 + OH. The H 2 + HCOH channel has 
a barrier, but it is below CH 3 + OH. The 1 CH 2 + H 2 0 channel has no barrier and 
is computed to be 1.4 kcal/mol below CH3 + OH. The CH3O + H and CH 2 OH + 
H channels are endoergic with respect to CH 3 + OH, but no intermediate barriers 
are expected. Thus, from this work the 1 CH 2 + H 2 0 , HCOH + H 2 , and CH 2 0 + 
H 2 channels are all found to be accessible from CH 3 + OH. 

Fig. 2 shows the energy as a function of rco for the addition of 1 CH 2 to H 2 0 , 
while the computed energies are given in Table IV. Two stationary points have been 
located along this minimum energy path. At long rco the dominant interaction 
is dipole-dipole, which results in an orientation with the dipole moments of the 
approaching molecules parallel to each other but pointing in opposite directions, 
as illustrated in Fig. 2 . At shorter rco there is a minimum on the PES with 
C 3 symmetry, followed by a saddle point with no symmetry for insertion of 1 CH 2 
into an OH bond of water. In order to characterize the minimum energy path, a 
gradient calculation was carried out starting at the saddle point and proceeding 
toward the minimum. CCI calculations were then carried out at the geometries 
corresponding to each step on the walk. This calculation defines the minimum 
energy path between the saddle point and minimum. It is not possible to do the 
same calculation for the portion of the PES connecting the minimum and the 1 CH 2 


8 



4- H 2 0 asymptote, since there is no well-defined way to follow the gradient uphill. 
In order to characterize this portion of the PES, calculations were carried out with 
the CH 2 and H 2 0 molecules fixed at their equilibrium geometries and oriented with 
the planes of the molecule parallel to each other and the CO bond perpendicular 
to both molecular planes in the orientation shown in Fig. 2. For this geometric 
orientation, rco was varied and the resulting energies at the CCI level are also 
included in Fig. 2. 

The main features of Fig. 2 are a shallow minimum followed by a small barrier 
to formation of CH 3 OH. However, the barrier is below the J CH 2 + H 2 0 asymp- 
tote and therefore the bottleneck on the vibrationally adiabatic curve is expected 
to occur in the entrance channel region. The main feature responsible for the en- 
trance channel bottleneck is the building in of bending modes, which arise from 
electrostatic (dipole-dipole, dipole-quadrapole, and quadrapole-quadrapole) inter- 
actions. In order to define this interaction, dipole and quadrapole moments (about 
the center of mass) were computed for 1 CH 2 and H 2 0 and are given in Table V. The 
experimental values [19] for H 2 0 (in a.u.) are /x = 0.73, Q xx — 1-955, Q yy = -1.859, 
and Q zz = -0.097; which are in reasonable agreement with the computed values. 

V. Conclusions. 

The potential energy surface (PES) for the CH 3 OH system has been charac- 
terized for the 1 CH 2 -f- H 2 0, H 2 + HCOH, and H 2 + H 2 CO product channels 
using complete-active-space self-consistent field (CASSCF) gradient calculations to 
determine the stationary point geometries and frequencies followed by CASSCF/ 
internally contracted configuration-interaction (CCI) calculations to refine the en- 
ergetics. 

The H 2 + H 2 CO, and H 2 + HCOH pathways have barriers located at 1.7 kcal/mol 
and - 5.2 kcal/mol with respect to CH3 -f OH. The a CH 2 + H 2 0 channel is found to 
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have no barrier in the absence of vibrational zero-point effects. However, the long 
range-interaction is dominated by a dipole-dipole term and the zero-point effects 
due to this interaction are expected to lead to a bottleneck on the vibrationally 
adiabatic minimum energy path. The a CH 2 + H 2 0 asymptote is computed to be 
1.4 kcal/mol below CH3 + OH. Thus, all three of these channels are expected to be 
accessible at moderate temperatures. 

From comparison of the computed energetics of the reactants and products to 
known thermochemistry it is estimated that the computed PES is accurate to 
± 2 kcal/mol. 

Acknowledgement. SPW was supported by NASA cooperative agreement number 
NCC2-478. 
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Table I. Computed Energetics for CHsOH 0 ’ 6 



calc. 

exp(OK) 

Harding et al. 

CH 3 + OH 


90.2 


HCOH-H 2 

85.0 


91.0 

HCOH + H 2 

71.1 


71.1 

a CH 2 + H 2 0 

88.8 

89. 5 C 

94.9 

ch 2 o-h 2 

91.9 


96.5 

ch 2 o + h 2 

18.1 

18.7 

16.4 

ch 3 oh 

0.0 

0.0 

0.0 

“ CASSCF/CCI with 

a [4s3p2dlf/3s2pld] basis set. 


b Relative energies in 

kcal/mol (including zero-point 

energy). 

c Using 3 CH 2 -> J CH 

2 separation of 9.0 kcal/mol. 
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Table Ila. Computed energies and zero-point corrections. 


CH 3 OH -+ CH 2 0 + H 2 



Energy 0 

zero-point energy 6 

A E' 

CH 3 OH 

-1 15. 51397(-. 54566) 

0.05258 

0.0 

CH 2 0-H 2 

-115. 35603(-. 39071) 

0.04407 

91.9 

ch 2 o+h 2 

-1 15. 47215(-. 50189) 

0.03768 

18.1 


a Energy in Eh- The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh- 
6 zero-point energy in Eh- 

c relative energy in kcal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table lib. Computed energies and zero-point corrections. 



CH 30 H 

CHi + H 2 0 



Energy“ 

zero-point energy 6 

A E c 

CH 3 OH 

-115.51397(-. 54566) 

0.05258 

0.0 

ch 2 -h 2 o 

-115. 36931 (-.40561) 

0.04366 

82.3 

ch 2 +h 2 o 

-115.35989(-. 39040) 

0.03884 

88.8 


a Energy in Eh- The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh- 
b zero-point energy in Eh- 

c relative energy in kcal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table lie. Computed energies and zero-point corrections. 



CH 30 H -> 

HCOH + H 2 



Energy 0, 

zero-point energy^ 

A E c 

CH 3 OH 

-1 15. 50589(-. 54231) 

0.05258 

0.0 

hcoh-h 2 

-1 15. 35803(-. 39640) 

0.04225 

85.0 

HCOH+H 2 

-115. 37823(-. 41465) 

0.03828 

71.1 


a Energy in Eh- The first energy is the CCI energy, while the energy in parenthesis 
includes a multi-reference Davidson correction and is with respect to -115. Eh* 
b zero-point energy in Eh- 

c relative energy in keal/mol including zero-point energy and a multi-reference 
Davidson’s correction. 
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Table III. Computed saddle point frequencies and rotational constants(cm *). 



CH 2 0-H 2 

HCOH-H 2 

ch 2 .h 2 o 

CH 2 -H 

UJ 1 

3195 

4130 

3860 

3723 

W 2 

2295 

3199 

3744 

3319 

W 3 

1740 

2323 

3187 

3222 

( 

1574 

1564 

3113 

2091 

^5 

1429 

1498 

1691 

1560 

U?Q 

1369 

1347 

1515 

1524 

CJ7 

916 

1291 

802 

1119 

U>8 

2877z 

1122 

610 

949 

0>g 

3278 

934 

337 

703 

^10 

1273 

622 

334 

438 


1211 

513 

161 

387 

^12 

1065 

1414i 

26 

1850z 

A 

3.345 

3.042 

4.117 

4.640 

B 

0.944 

0.839 

0.358 

0.496 

C 

0.863 

0.764 

0.349 

0.478 
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Table IV. Computed energies and zero-point corrections 0 . 




CH 2 + H 2 0 ch 2 

-H 2 0 


r co 

CAS Energy 

Cl Energy 6 

zero-point energy 

total 0 

3.834 

-115.00924 

-115. 36931(-. 40561) 

0.04366 

(-.36225) 

3.929 

-115.01305 

-115.36889(-.40409) 

0.04449 

(-.35960) 

4.076 

-115.01765 

-115.36968(-. 40420) 

0.04550 

(-.35870) 

4.207 

-115.02142 

-115. 37055(-. 40452) 

0.04546 

(-.35906) 

4.310 

-115.02426 

-115.37130(-. 40480) 

0.04523 

(-.35957) 

4.386 

-115.02635 

- 1 15. 37189(-. 40503) 

0.04500 

(-.36003) 

4.437 

-115.02785 

-115.37230(-.40516) 

0.04482 

(-.36034) 

4.471 

-115.02890 

-115. 37254(-. 40517) 

0.04470 

(-.36047) 

4.528 

-115.03024 

-1 15. 37237(-. 40450) 

0.04421 

(-.36029) 

5.0 

-115.02970 

-11 5. 36971 (-.401 36) 



5.5 

-115.02892 

-115.36727(-. 39843) 



6.0 

-115.02797 

-115.36541(-. 39630) 



6.5 

-115.02710 

-115.36400(-.39475) 



7.0 

-115.02635 

-1 15. 36294(-. 39358) 



8.0 

-115.02526 

-115. 361 51 (-.39206) 



9.0 

-115.02464 

-115.36076(-.39127) 



10.0 

-115.02433 

-115.36040(-. 39091) 




a All energies in Eh* *co is in a. u. 
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b The first energy is the CCI energy, while the energy in parenthesis includes a 
multi-reference Davidson correction and is with respect to -115. Eh- 
c Energy with respect to -115. Eh* The energy includes the zero-point correction 
and a multi-reference Davidson correction. 
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Table V. Computed dipole and quadrapole moments a . 


fl Qxx Q yy Qxz 

HoO 0.768 -1.987 1.955 0.032 

CH 2 0.679 1.187 0.496 -1.683 


a Properties are in a.u. Quadrapole moment is with respect to the center of mass. 
The molecule is in the YZ plane with the C 2 axis in the Z direction. 
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Figure Captions. 


Fig. 1. Schematic diagram of the potential energy surface for CH3 + OH. The 
location of the CH3 + OH asymptote with respect to CH3OH is taken from experi- 
ment, while the locations of the CH 2 OH + H aand CH3O + H asymptotes are from 
previous calculations. 

Fig. 2. The potential for 1 CH 2 + H 2 0 from CCI calculations along the CASSCF 
minimum energy path. See the text. 
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Computed Barrier Heights for 





- p. ? 

N 93'- 17409 


H + CH 2 0 <-► CH s O «-* CH 2 OH 


Stephen P. Walch a 
ELORET Institute 
Palo Alto, Ca. 94303 

Abstract. The barrier heights (including zero-point effects) for H + ChEO — * CH 3 O 
and CH 3 O —+ CHoOH have been computed using complete active space self consis- 
tent field (CASSCF)/ gradient calculations to define the stationary point geome- 
tries and harmonic frequencies and internally contracted configuration-interaction 
(CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 
30.1 kcal/mol, respectively. The former barrier height compares favorably to an 
experimental activation energy of 5.2 kcal/mol. 
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The CH 3 O radical is important in combustion and atmospheric chemistry [ 1 ]. In 
addition, recent experiments have studied individual rovibrational states of CH3O 
above the H-CH 2 0 dissociation threshold by stimulated emission pumping [3]. The 
interpretation of these experiments depends critically on the barriers to decom- 
position of CH 3 O to H + CH 2 0 and rearrangment to CH 2 OH. Previous studies 
of these barrier heights include the work of Saeb 0 , Radom, and Schaefer [2], who 
used Mpller-Plesset perturbation theory through third order with small basis sets 
up through 6-31G**. More recently, Page, Lin, He, and Choudhury [4] reported 
somewhat more accurate calculations using multireference configuration-interaction 
(MRCI) with polarized triple zeta basis sets. In the present paper more accurate 
calculations are reported which make use of larger basis sets and extensive Cl. 

Two different basis sets were used in this work. For the complete active space 
self consistent field (CASSCF) gradient calculations the polarized double zeta set 
of Dunning and Hay [ 5 ] was used. The basis set for C and 0 is a (9s5p)/[3s2pj 
basis augmented by a single set of 3 d functions with exponents of 0.75 and 0.85 
for C and 0, respectively. The H basis is (4s) / [ 2 s] augmented with a single set of 
2 p functions with exponent 1 . 00 . The basis set used in the Cl calculations is the 
Dunning correlation consistent triple zeta double polarization basis set [ 6 ]. This 
basis is [ 4 s 3 p 2 dlf] for C and 0 and [ 3 s 2 pld] for H and is described in detail in Ref. 
6 . 

The calculations were carried out in C* symmetry for a wavefunction of 2 A 1 
symmetry. From the Hessian matrix in the diagonal representation, the station- 
ary points obtained in C 5 symmetry are the minimum for CH 2 0 + H, the saddle 
point for CH 2 0 + H CH 3 O, the CH 3 O minimum, and the saddle point for 
CH 3 O — > CH 2 OH. For CH 3 O the 2 A' state corresponds to one component of the 
2 E state in C 3 V symmetry. According to the calculations of Saeb 0 et al. [ 2 ], the 
Jahn-Teller stabilization energy in this system is 0.56 kcal/mol and the Jahn-Teller 
splitting is 0.12 kcal/mol, with the 2 A' component lower. CH 2 OH is found to have 
no symmetry and the use of a mirror plane of symmetry would not be appropriate. 
However, the separation between CH 3 O and CH 2 OH has been accurately computed 
by Bauschlicher, Langhoff, and Walch [7] and further calculations for CH 2 OH were 
not carried out here. 
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The CASSCF calculations consisted of five electrons and five orbitals. The active 
electrons for CH 3 0 included the CO bond pair, one CH bond pair (the CH bond in 
the molecular plane), and the a' 0 2p like orbital. These orbitals correspond to the 
CO o and 7 r bonds plus a H Is orbital for CH 2 0 + H. The remaining electrons, which 
are inactive, include the other two CH bond pairs, the a" 0 2p like orbital and the 
0 Is, 0 2s, and C Is like orbitals. In generating the set of reference configurations 
for the subsequent internally contracted Cl (CCI) calculations, no more than two 
electrons were permitted in the weakly occupied CASSCF orbitals. All but the 0 Is 
and C Is electrons were correlated in the CCI calculation. A multi-reference analog 
of Davidson’s correction [8] was added to the CCI energies. 

The CASSCF /gradient calculations used the SIRIUS/ABACUS system of pro- 
grams [9], while the CCI calculations were carried out with MOLPRO [10,11]. The 
calculations were carried out on the NASA Ames Cray Y-MP. 

Table I gives the harmonic vibrational frequencies, rotational constants, and ge- 
ometric parameters for the stationary points considered in this work. Here H a is 
the H of the spectator CH bonds, while H& is the hydrogen which transfers during 
the reaction. The geometries and frequencies in Table I are quite similar to those 
obtained by Page et al. [4]. This is expected since a CASSCF wavefunction was 
used in Ref. 4 . However, in Ref. 4 the CO o bond was not correlated. For CH 3 0 
this has the expected effect that the CO bond is elongated by ^ 0.03 A compared to 
the calculations in Ref. 4 . At the H-CHoO saddle point the CO bond is elongated 
by ^ 0.02 A and the C-H b bond is shorter by 0.05 A as compared to Ref 4. 
There are larger differences between this work and the results of Saebp et al. Here 
the geometry for CH 3 O is similar except for the expected elongation of the C-H& 
and CO bonds. However, there are significant differences in the geometry of the 
H-CH 2 O saddle point. In particular, the C-H& bond length is shorter by ~ 0.18 A 
. This result is consistent with the expectation that saddle point geometries will be 
more sensitive to electron correlation effects than will stationary points. 

Table II gives the computed energies and zero-point corrections. The energies 
are from the CCI calculations, while the zero-point corrections are derived from the 
CASSCF harmonic frequencies. From Table II it is seen that the relative energies are 
only very slightly changed 0.1 kcal/mol) in going from the [ 4 s 3 p 2 dlf/ 3 s 2 pld] 
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basis set to the [5s4p3d2f/4s3p2d] basis set. Table III compares the computed 
barrier heights obtained in the present calculations to those obtained by Saebp et 
al. [2] and by Page et al. [4]. Here it is seen that the barriers obtained by Saebp 
et al. are too large by « 6 kcal/mol., while the barrier for H + CH2O — > CH3O 
obtained by Page et al. is 2.3 kcal/mol larger than in the present calculations. Our 
computed barrier height is 5.6 kcal/mol for H + CH2O — > CH3O and 5.2 kcal/mol 
for D + CH2O — > CDH2O. This result is in good agreement with experimental 
estimated activation energies of 5.2 kcal/mol for H + CH2O and 3.9 kcal/mol for 
D + H 2 CO [12]. 

It is concluded that the present calculations give energetics for the CH3O sys- 
tem which are of chemical accuracy ( in this case the estimated error bars are 
± 1 kcal/mol). This is in contrast to results obtained by Saebp et al. using Mpller- 
Plesset perturbation theory, where the errors in the computed barrier heights are 

6 kcal/mol. 

Acknowledgement. SPW was supported by NASA cooperative agreement NCC2- 
478. 
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Table I. Computed stationary point harmonic frequencies (cm J ), rotational 
stants (cm -1 ), bond lengths (A), and bond angles (degrees). 



CH 3 0 

CH 3 0 - 

■+ CH 2 OH 

H-CH 2 O 

H + CH 2 0 


3216 

3261 


3190 

3174 

0J2 

2921 

2378 


1657 

1766 


1608 

1556 


1471 

1589 

(x/4 

1471 

1101 


1220 

1213 

^5 

1103 

1037 


605 


u; 6 

997 

2254 i 


1435 i 


U>7 

3284 

3391 


3297 

3272 

Ct>8 

1491 

1159 


1281 

1304 

U? 9 

1050 

1015 


770 


A 

5.180 

5.718 


3.927 

9.548 

B 

0.885 

0.929 


1.007 

1.263 

C 

0.884 

0.900 


0.962 

1.115 

*c-H a 

1.092 

1.086 


1.094 

1.097 

*C-H b 

1.122 

1.275 


1.668 


O 

1 

O 

1.412 

1.418 


1.259 

1.226 

ro-H t 


1.237 


2.292 


Z H a CH a 

111.5 

119.1 


117.1 

117.2 

Z H a CH 6 

108.2 

117.1 


91.2 


Z H a CO 

111.7 

116.5 


120.3 

121.4 

Z H 6 CO 

105.3 

54.4 


102.2 



con- 
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Table II. Computed energies and zero-point corrections. 
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Table III. Comparison of computed barrier heights (kcal/mol). 




Saebp et al. a 

Page et al. fr 

PW C 

Exp. 

H + CH 2 0 

— CH 3 0 

12.5 

8.0 

5.6 

5.2 d ,6.2' 

D + CH 2 0 

-+ cdh 2 o 



5.2 

3.9 d ,4.8 


CH 3 0 -> CH 2 OH 36.0 30.1 

a Reference 2. 
b Reference 4. 
c Present work. 

d Estimated activation energy from a room temperature rate constant for D + 
H 2 CO, an assumed preexponential factor of 2 x 10 13 cm 3 mol 1 sec 1 , and a zero- 
point energy difference between D T H 2 CO and H 4- H 2 CO of ~ 1.3 kcal/mol. Ref. 
12 . 

e Estimated threshold energy from an RRKM analysis. Ref. 12. 
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